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ABSTRACT 

Context. Emission and absorption line profiles which are formed in the interaction of pre-main sequence (PMS) stars and their 
circumstellar environment are found to be variable at various timescales. 

Aims. We investigate the patterns and timescales of temporal line profile variability in order to explore the dynamical circumstellar 
environment of the PMS Herbig Ae star V35 1 Ori. 

Methods. We obtained 45 high-resolution (R~28,000) spectra of V351 Ori at timescales of hours, days, and months. We analysed the 
Ha line profiles and also examined the Hp, NaDl and NaD2 line profiles to explore the nature of the spectroscopic variability. 
Results. The Ha line profiles showed strong variations over all timescales. The shape of the profiles changed over timescales of a 
day. Single as well as simultaneous event(s) of blue-shifted and red-shifted transient absorption components (TACs), i.e. signatures of 
outflow and infall, were also observed in the Ha line profiles. The shortest period of variation in the TACs was < 1 hour. All transient 
absorption events were found to decelerate with a rate of a few to fractions of m s~ 2 . The depth and width of the TACs were also 
changing with time. The presence of elongated red-shifted components at some epochs supports the episodic nature of accretion. 
Conclusions. Variable emission and absorption components detected in Ha line profiles show the dynamic nature of interaction 
between V35 1 Ori and its circumstellar environment. The Ha non-photospheric profiles of the star most probably originate in the disk 
wind. Episodic accretion of gaseous material at a slow rate and outflow of clumpy gaseous material are still occurring in V351 Ori 
at an age of ~6.5 Myr. Dynamic magnetospheric accretion and disk wind emerge as the most satisfactory model for interpreting the 
observed line profile variations of V35 1 Ori. 
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1. Introduction 

Herbig Ae/Be (HAeBe) stars are revealed pre-main sequence 
(PMS) emission line stars with masses of ~1.5 to 10 M Q and 
typical ages of ~10 6 years with spectral types earlier than F0. 
Balmer emission lines of hydrogen and infrared (IR) excess 
emission produced by thermal emission of circumstellar dust 
grains are the prominent characteristics of HAeBe stars. Both 
of these signatures originate in the circumstellar material around 
these stars and the interaction between them. It has been well 
established that circumstellar material, i.e. dust and gas around 
the HAeBe stars, are distributed in disks (e.g. Maheswar et al. 
(2002) and references therein). Significant variations in shape 
and intensity of permitted and forbidden emission line profiles, 
originating in the accretion or outflow events of the interaction 
between the central star and its surrounding material, are ob- 
served in many HAeBe stars ( |Hamann & Persson|1992 ; Hamann 
[1994 



Reipurth et al.| T996J. Optical spectroscopic studies of 



emission line profiles of a few HAeBe stars are reported in the 
literature in e.g. |Grinin et al.| (|1994) ; Ide Winter et al.| d l999 



Grinin et al.|(|2001);|Natta et al. (2000b;|Mora et al.|d2002| 12004 



|Guimaraes et al.| ( |2006) . Natta et aL| ( |2000[ l reported the evidence 
for episodic rather than continuous accretion in Herbig Ae star 
(2002, 2004) reported transient absorption 
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components (TACs) in several Balmer and metallic lines of a few 
Herbig Ae type stars, which are supposed to be created by mov- 
ing gaseous blobs around these stars. These studies discussed 
the dynamical circumstellar environment around HAeBe stars 
which show strong Ha emission. 



Manoj et al. ( 2006[ ) showed that the emission line activity of 
HAeBe stars substantially decreases by ~3 Myr, which agrees 
well with the inner disk survival timescale i.e ~3 Myr as sug- 



gested by Hernandez et al. (2005 2009) for the early spectral 
type stars such as B, A, and F types. The dissipation of optically 
thick disks beyond ~3 Myr to the onset of young main sequence 
stars with a debris disk can be caused by various disk disper- 
sal processes such as photoevaporation, planet formation, etc. 
HAeBe stars with weak emission lines and moderate IR excess 
can be considered as the precursor of the young main sequence 
stars with IR excess and no gaseous emission lines, e.g. (3 Pic 
and Vega. Understanding the nature of accretion and circumstel- 
lar material of comparatively less active HAeBe stars of age > 
3 Myr are important because these processes eventually set the 
stage for planet formation. However, detailed studies of weak 
emission line HAeBe stars that cover hours, days, and months 
are not available in the literature. 

V351 Orionis (V351 Ori, HIP 27059, HD 38238) is a weak 
Ha emission line Herbig Ae star that showed photometric and 
spectroscopic variabilities over various timescales in e.g. |van 
|den Ancker et al.| (jl996); Balon a~et~aLl d2002| ) and references 
therein. The classification of V35 1 Ori as a Herbig Ae star has 
been confirmed by several authors, e.g van den Ancker et al. 
" Vieira et al.| ( |2003[ ); |Hernandez et al.|d2005| ). |Grady et al. 



(1996 



( 1996) analysed the ultraviolet (UV) spectra of V351 Ori and did 
not find any absorption component, i.e. signatures of stellar wind 
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or accretion, in the spectra, van den Ancker et al. ( 1996) pre- 
sented photometric and spectroscopic observations of V35 1 Ori 
and found that it transformed into an almost non-variable star 
from a strong photometric variable within a short period of ~ 14 
years. |van den Ancker et al.| reported that the Ha profile of 
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V351 Ori was of an inverse P Cygni type and also showed a high 
infrared excess emission. Inverse P Cygni profiles are signatures 
of the accretion of circumstellar material onto the central object 
and are characterised by blue-shifted emission with red-shifted 
absorptions or systematically enhanced red shifted absorptions. 
V351 Ori shows 5 Scuti type pulsations (Marconi et al. 2000), 
which are another characteristic of PMS stars of mass > 1.5M 
as they cross the pulsation instability towards their contraction 
to the main sequence. Balon a~et al.| < |2002] > obtained simultane- 
ous multi-colour photometric and spectroscopic observations of 
V351 Ori to investigate the nature of the pulsation of V351 Ori. 
Balo na et al.| also reported a variable inverse P Cygni type Ha 
profile with a variation on a timescale of a day with complete ab- 
sence of emission at some epochs. Vieira et al. (2003) observed 
an inverse P Cygni type Ha emission line profile in V351 Ori 
and also suggested that V35 1 Ori might be an evolved Herbig Ae 
star, because it shows weak emission lines. Herna ndez et al.| 
(2005 ) also reported a low equivalent width of Ha i.e EW^[Ha ] 
=-0.9 A and associated V351 Ori with Ori OBlbc region. These 
signatures hint at an active interaction of the weak Ha emission 
line young star V35 1 Ori with its circumstellar environment. 

In this work we present high-resolution spectroscopic ob- 
servations of V351 Ori along with optical BVRI observations. 
We discuss the dynamic circumstellar gaseous environment of 
the star in the light of available models of star-disk interac- 
tion. Observations and data reduction are discussed in Sect. [2] 
In Sect. [3] we present the details of the spectra and the analysis 
of the TACs detected in the spectra. A discussion of the kine- 
matics of the observed TACs, an estimate of the disk mass and 
a qualitative interpretation of the Ha line profile variation are 
presented in Sect. |4] We present our conclusion in Sect. [5] 

2. Observations and data reduction 

Fortyfive high-resolution spectra of V351 Ori were obtained dur- 
ing October 2008 to April 2009 with a fibre-fed cross-dispersed 
echelle spectrometer (R ao et al.|200 5 ) at the 2.3 m Vainu Bappu 
Telescope of Vainu Bappu Observatory (VBO), Kavalur, IndiaH 
The log of the observations is given in Table [T] The resolving 
power of the spectrometer set-up used for observation is ~28,000 
and the attached 2Kx4K CCD gives a dispersion of ~0.025 A per 
pixel. The typical integration time of each spectrum was 45 min- 
utes. Bias subtraction, flat field correction and scatter light re- 
moval were carried out for all the spectra with the standard tasks 
available in the Image Reduction and Analysis Facility (IRAF^ 
The Th-Ar lamp was used for wavelength calibration. The fina 
spectra that were used for our analysis cover the wavelength re- 
gion from A4250 A to /I7900 A. The CCD used for observation 
did not cover the full echelogram, and as a consequence there are 
some gaps in the wavelength coverage of the observed spectra. 
However, we optimised the CCD position accordingly to cover 
the wavelength regions of Ha, H/3, NaDl, NaDl and several 
other photospheric lines. Atmospheric lines were also used to 
check the wavelength calibration. These wavelength calibrations 
of the final spectra are accurate up to +0.01 A. 

V351 Ori was observed on 11 November 2009 with the 
Bessell broad band filters B (5x5 s), V (4x2 s), R (4x1 s) and 
/ ( 5x1 s) of the Himalayan Faint Object Spectrograph Camera 
(HFOSC) mounted on the 2m Himalayan Chandra Telescope 



Table 1. Observation log of the echelle spectroscopy of 
V351 Ori 



1 http://www.iiap.res.in/vbo_vbt 

2 The IRAF software is distributed by 



Date 


TT"\ T A C A AAA 

JD - 2454000 




Ha profile type" 


2008-10-27 


767.28-767.41 


2700 x 4 


IV-R,inverse P Cygni 


2008-10-30 


770.38 


2700 x 1 


II-R, double peak 


2008-12-11 


812.32 


2700 x 1 


IV-R 


2008-12-28 


829.10-829.32 


2700 x 6 


IV-R 


2008-12-29 


830.10-830.37 


2700 x 8 


IV-R 


2009-01-17 


849.11 -849.22 


2700 x 4 


III-R, double peak 


2009-01-18 


850.10-850.32 


2700 x 4 


III-R 


2009-01-19 


851.09-851.35 


2700 x 5 


III-R 


2009-01-20 


852.07-852.18 


2700 x 2 


III-R 


2009-03-02 


893.11 -893.15 


2700 x 2 


III-R 


2009-03-03 


894.20 


2700 x 1 


II-R 


2009-03-29 


920.16 


2700 x 1 


IV-R 


2009-03-30 


921.07-921.11 


2700 x 2 


II-R 


2009-03-31 


922.08 


2700 x 1 


III-R 


2009-04-26 


948.07 


2700 x 1 


III-R 


2009-04-27 


949.07 


2700 x 1 


II-R 


2009-04-28 


950.07 


2700 x 1 


II-R 



a. Based on the classification scheme adopted from Reipurth et al. 1 1996). Type 
II-R profile represents a double peak profile where the red-sided peak exceeds 
the strength of the blue-sided peak by half. Type III-R also represents double 
peak profile with the red-sided peak less than half the strength of the blue-sided 
peak. Type IV-R represents the inverse P Cygni profile. 



(HCT) of the Indian Astronomical Observatory (IAO), Hanle, 
IndiqJ Data were reduced with the standard tasks in IRAF. 



3. Results and analysis 

3.1. Description of echelle spectra 

HAeBe stars show both short and long term variabilities in their 
spectra. To explore the short to long term variations of the line 
profiles we obtained high-resolution spectra on timescales as 
short as an hour to the longest coverage of ~7 months. All 45 
spectra cover the wavelength regions of Ha, HQ, NaDl, NaDl 
and several other photospheric lines. Typical signal- to-noise ra- 
tios (S/N) are < 10 to 12 at HB, 15 to 30 at NaDl and NaDl and 
> 30 to 60 at Ha. We note that [O I] ,16300 A emission line was 



present in several spectra. IBalona et al. ( 2002 1 also reported the 



is 

Astronomy Observatory under contact 
Foundation, http : //iraf .noao . edu/ 



the National Optical 
with National Science 



detection of a narrow emission line of [O I] /I6300 A at a veloc- 
ity of +7.5 km s -1 . However, the line centre of [O I] /16300 A 
in our observed spectra remained constant at laboratory wave- 
length, as did the other telluric lines throughout our observing 
runs. We therefore conclude that the [O I] /16300 A emission line 
in our observed spectra is of telluric origin. We did not find any 
detectable emission in the HB, NaDl and NaDl line profiles, 
including the other photospheric lines typical to an A7 type stars 
within the noise level of our observed spectra. Though the Ha 
line profiles were dominated by variable absorption components, 
but nominal to significant variable emission components were 
also observed. Below we describe in detail the behaviour of dif- 
ferent line profiles. To compare the spectra obtained in different 
epochs one needs to normalise the spectra in a similar fashion. 
The normalisation of the echelle spectra by fitting the continuum 
is tricky, because each order displays a limited wavelength cov- 
erage (in this case ~40 A). However, with the synthetic spectra 
of the star as the guiding standard, a better fit to the continuum 
and accordingly a common normalisation of all spectra can be 
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obtained. As we shall discuss below, synthetic spectra are also 
essential to identify the circumstellar component in the line pro- 
files that are dominated by absorption components. 

3.2. Stellar parameters and synthetic spectra 



Balona et al. (2002) estimated the effective temperature (T e //), 
the rotational velocity (v sin i), and the log g (logarithm of 
the surface gravity) of V351 Ori in order to generate the syn- 
thetic spectra of the star. Ripepi et al. (2003) estimated effective 
temperature (T c yy ~7425-7600 K) of the star using the pulsa- 
tion models. The average effective temperature of Ripepi et al 



matches quite well with the estimated value of 7500 K by |Balona 
|et al.| Recent studies on the stellar parameters of HAeBe stars 
showed that the typical value of log g is ~4.0, which also gives 
a reasonable fit to the photospheric wings of the Ha and H/3 
profiles of emission line HAeBe stars. We used the model at- 
mosphere from the Kurucz model atmosphere databas^jfor an 
A7 II I star with T e ff =7500 K as adopted from Balo na et al 
d2002|l and log g= 4.00. We used v sim'=102 km s~ 1 ( |BaIona 
|et al.|[2002| and a micro turbulent velocity of 2 km s~ ! as sug- 
gested by Dunkin et al. ( 1997 1. Detailed synthetic spectra were 
computed with the interactive data language (IDL) interface 
SYNPLOT (I. Hubeny, private communication) to the spectrum 
synthesis programme SYNSPEC |Hubeny & Lanz| ( |2000) l, utilis- 
ing the adopted Kurucz model atmosphere. The linelist provided 
in SYNSPEC is used to generate the synthetic spectra. 

We convolved the synthetic spectra with the rotational pro- 
file of V35 1 Ori and the instrumental profile of the echelle spec- 
trometer with SYNPLOT and ROTINS0 In Figs, [I] and [2] we 
display some of the representative photospheric line profiles and 
overplot the convolved synthetic spectra. Because the observed 
photospheric and synthetic line profiles show a good match over 
long wavelength ranges, we did not attempt to fit the photo- 
spheric lines rigorously to obtain the stellar parameters. We used 
the convolved synthetic spectra to identify unblended photo- 
spheric lines to obtain the radial velocity of the the star. We used 
Fe I, Fe II, and Ca I lines in a wavelength range from ,15000 
to /16500 A [such as Fe II ,14923.827 A, Fe I /15001.862 A, 
Fe I ,15429.696 A, Fe II ,15534.847 A, Fe I ,15615.644 A, 
Ca I ,16102.723 A, Ca I ,16122.217 A, Ca I ,16439.075 A] to 
measure the radial velocity of the star. We estimate the average 
heliocentric radial velocity of the star as + 11 +3 km s _1 , which 
is similar to the estimated radial velocity i.e. + 13 km s _1 by 



Balona et al.| We synthesised the photospheric Ha line profile 
following the same procedures. 

3.3. Ha line profiles 

We display all the Ha emission line profiles in the rest frame ve- 
locity of V35 1 Ori and overplot the photospheric components in 
Fig|3ja) to Fig[6|a). We also mark the spectra by their date of ob- 
servation and the modified heliocentric Julian day (MHJD) de- 
fined by 2454000-heliocentric Julian day (HJD). We subtract the 
photospheric component from the observed spectrum and show 
the residual spectrum in Fig[3jb) to Fig [6jb). The zero level in 
the residual spectrum is described as the zero absorption level. 
As discussed before, profile variations on a timescale of days to 
months are quite common to V351 Ori. We also observed both 
blue- and red-shifted absorption components of various widths 



4 http://kurucz.harvard.edu/ 
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Fig. 1. Average photospheric line profiles of V351 Ori on 
29 December 2009 (dotted green line). Synthetic spectra in 
the same wavelength range are overplotted (solid blue line). 
The prominent photospheric lines are Fe I 5429.696, 5434.523, 
5445.042, 5446.871, 5446.916, 5455.441 and 5455.609 in the 
first row, Sc II 5526.790, Mg I 5528.405 and Fe II 5534.847 
in the second row, and Fe I 6430.844, Fe II 6432.680, Ca I 
6439.075, and Fe II 6456.383 in the third row. (This figure is 
available in colour in electronic form.) 



and depths in several epochs. These features are also seen in 
other HAeBe stars su ch as UX Or i , BF Ori, SV Cep, WW Vul, 
and XY Per, 



see e.g. 



Mora et al. (20021^004) and references 



therein. These features are called transient absorption compo- 
nent (TAC) because they are sporadic in nature and last for a few 
days only. Transient absorption components can be divided into 
two categories, i.e blue-shifted absorption components (BACs) 
and red shifted absorption components (RACs). The chemical 
compositions and kinematics of the TACs have been discussed 
by se veral studies, e.g. |Natta et al] ( |2000| i; |Mora et al] ( |2002 



2004). We discuss the profile variations in detail in the subse 



quent sections. We quote the velocity of the broad wings from 
the residual spectra and that of the TACs from the observed nor- 
malised spectra. 

3.4. Line profile variations 



Previous stu dies of e.g. |van d en Anc ker et al.| ( |1996| ), |Balona| 
et al. ( 2002} , and |Vieira et~aL ( |2003 1 reported that the average 
Ha emission line profiles of V351 Ori were inverse P Cygni 
types. From Table. [T]it is evident that V351 Ori showed various 
profile shapes such as II-R and III-R, including inverse P Cygni 
types during our observation. Nightly variations of profiles from 
an average behaviour at a timescale of a day were also ob- 
served in several epochs. The type II-R profile of 2 March 2009 
evolved into a type III-R on 3 March 2009. Another incident 
of a rapid variation in profile shapes was observed during 29 
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Fig. 2. Average photospheric line profiles of V351 Ori of 17 to 
19 January 2009 {dotted green line). Synthetic spectra in the 
same wavelength range are overplotted (solid blue line). (This 
figure is available in colour in electronic form.) 
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Fig. 3. (a), //a line profiles {green solid lines) of V351 Ori on 
27th October 2008. The synthetic Ha line profiles (blue dashed 
lines) are overplotted in each spectrum. Modified Julian days 
are displayed against the respective spectrum. The vertical red 
dashed line represents the rest frame velocity of the star. (b). 
Residual spectra of Ha line profiles on same epochs. Horizontal 
blue dashed lines represent the zero absorption levels. (This fig- 
ure is available in colour in electronic form.) 



to 31 March 2009. Three different types of Ha profiles, namely 
IV-R, II-R, and III-R, were observed on these consecutive nights. 
The flip-over from type II to type III profiles, i.e the change in 
strength of blue- and red-shifted emission peaks in a double- 
peak profile, was also reported by Vieira et aL] ( 2003 i in PDS 018 
and PDS 024, but on longer timescales that were separated 
by months to year respectively. We also observed two promi- 
nent incidents of infall and outflow in the Ha profiles on 28 to 
29 December, 2008 and 17 to 20 January, 2008 respectively. The 
infall event was seen in the H/3 profile as well (Fig. [TT] ). We 
also checked a few photospheric metallic lines in these epochs 
that were separated by nearly two weeks (Figs. T]and 2]) and did 
not find any variation in those profiles, which also indicates that 
there was no significant variation in the stellar photosphere dur- 
ing the infall and outflow events. 

We calculated the mean and temporal variance profile of Ha 
line profiles to investigate the characteristics of the variations. 
We calculated the mean profile by taking the average normalised 
intensity in each velocity interval of a given profile over all ob- 
served profiles. We show the average normalised Ha line profile 
of V351 Ori in Fig. [7] As suggested by |Johns & BasrT| (1995) , 
we define the temporal variance profile at each velocity interval 
as 



Z 



n-l 



where n is the number of spectra, 7 lv is normalised intensity at 
a given velocity (v) in each spectrum (i), and I v is mean nor- 
malised intensity at a given velocity (v) over all observed spec- 
tra. The normalised variance profile (shown as shaded area in 



Fig. [7]) can be obtained by dividing the temporal variance profile 
by the average profile. The normalised variance profile indicates 
that the behaviour of an individual spectrum observed at differ- 
ent epochs deviates substantially from the average profile. Most 
of the variability in Ha profile is on red side with two clearly 
distinguished features. There are significant variations on blue 
side as well, but the profile shows a single peak. A couple of 
outflow events were observed at different epochs separated by 
days to months, but they seem to have an approximate average 
velocity of -100 km s _1 . We also observed simultaneous infall 
and outflow events of different durations at different velocities. 
It is obvious from Fig.|7]that the infall events occurred at two dis- 
tinctly different velocities of e.g. ~60 km s _1 and ~150 km s . 
In emission line stars, circumstellar emission first fills the pho- 
tospheric absorption components and the residual flux emerges 
as emission from the star. If there is strong emission, the reduc- 
tion of emission flux to cover up the photospheric absorption 
is not significant. However, for weak emission, where the Ha 
line profile is dominated by absorption rather than emission, a 
subtraction of photospheric component is necessary to achieve 
the true contribution caused by emission. Red-shifted narrow ab- 
sorption components in Fig.[3|a) to Fig. |6ja), which indicate in- 
fall of matter, are prominently visible beyond the photospheric 
profiles. Accordingly the subtraction of photospheric profiles are 
not crucial for comparative studies of these features. However, 
the broad blue- and red-shifted emission/absorption components 
at the wings sometimes appeared at the same level as the pho- 
tospheric profiles. To make sure that we were not looking at the 
photospheric, but at the circumstellar components, we subtracted 
the photospheric spectra from the observed normalised spectra 
and plotted the profiles in Fig.[3jb) to Fig.|6jb). The residual pro- 
files also provide a better estimate of maximum velocity reached 
at the wings of the profiles, which are not contaminated by pho- 
tospheric components. We define the maximum velocity of the 
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Fig. 4. (a). Ha line profiles {green solid lines) of V351 Ori of 
some days in December 2008. Synthetic Ha line profiles {blue 
dashed lines) are overplotted in each spectrum. Modified Julian 
days are displayed against the respective spectrum. The verti- 
cal red dashed line represents the rest frame velocity of the star, 
(b). Residual spectra of Ha line profiles in the same epochs. 
Horizontal blue dashed lines represent the zero absorption lev- 
els. (This figure is available in colour in electronic form.) 



Fig. 5. (a). Ha line profiles {green solid lines) of V351 Ori of 
some days in January 2009. Synthetic Ha line profiles {blue 
dashed lines) are overplotted in each spectrum. Modified Julian 
days are displayed against the respective spectrum. The verti- 
cal red dashed line represents the rest frame velocity of the star, 
(b). Residual spectra of Ha line profiles in the same epochs. 
Horizontal blue dashed lines represent the zero absorption lev- 
els. (This figure is available in colour in electronic form.) 
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Fig. 6. (a). Ha line profiles (green solid lines) of V351 Ori of 
some days in March and April 2009. Synthetic Ha line profiles 
{blue dashed lines) are overplotted in each spectrum. Modified 
Julian days are displayed against the respective spectrum. The 
vertical red dashed line represents the rest frame velocity of the 
star. (b). Residual spectra of Ha line profiles in the same epochs. 
Horizontal blue dashed lines represent the zero absorption lev- 
els. (This figure is available in colour in electronic form.) 



profile at wings as the full-width-at-zero intensity of the resid- 
ual spectra. We also calculate the average and variance profiles 
using residual spectra. The average profile in Fig [8]is a double- 
peak profile with a strong blue to red asymmetry, i.e. III-R type, 
which is different from its counter part in Fig [7] The red- shifted 
absorption near the centre is similar to that of Fi g [7] and there is 
no change in the variance profiles in Figs |7j and [8] In Fig. [HJ the 
two peaks of the red-shifted variance profiles do not coincide 
with red-shifted emission peak of average profile of the resid- 
ual spectra, which supports our assumption of emission in most 
epochs. Double peak profiles are supposed to form in stellar or 
disk winds contrary to inverse P Cygni profiles, which originate 
in the infall of materials. We shall discuss below the effects of 
accretion, winds, etc. on the profile shapes . Significant changes 
in the average profiles in Figs. [7] and [8] show that photospheric 
profile subtraction is necessary at least for weak Ha emission 
line stars. 
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Fig. 7. Average Ha line profile (solid line) and variance profile 
(grey shaded area) of V35 1 Ori obtained using all spectra. 
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Fig. 8. Average Ha line profile (solid line) and variance profile 
(grey shaded area) of V35 1 Ori obtained using all residual spec- 
tra. 



3.5. Profile variations on 27 and 30 October, 2008 

We observed a double-peak profile, e.g. a type III-R profile, 
on 27 October 2008 with moderate blue-shifted and weak red- 
shifted emission components just above the photospheric spectra 
(Fig.[3ja)). The blue-shifted emission component did not change 
much, but the red-shifted component gained strength and be- 
came comparable to its blue-shifted counterpart, i.e. there were 
nearly symmetric broad emission components on both sides of 
the central minima, and thus the profile shape evolved into a type 
II-R on 30 October. The maximum average velocity of the pro- 
file was 350 km s _1 and ~250 km s _1 on 27 October, and 

300 km s _1 and ~280 km s _1 on 30 October at the blue and 

red wings respectively. We observed simultaneous inflow and 
outflow events on 27 October. The central broad asymmetric ab- 
sorption feature (~40 km s~') on 27 October was extended more 
towards the red region. It consisted of two barely resolved com- 
ponents with centres at ~6 km s~'and ~50 km s _1 . This broad 
absorption might be a combination of two absorptions of nearly 
similar depths because of both the disk and the infall of the cir- 



R. Choudhury et. al.: Variable circumstellar activity of V351 Ori 



7 



cumstellar material. The width of the central absorption was re- 
duced on 30 October and seemed to represent a single absorption 
component at ~9 km s _1 , which was nearly symmetric with the 
rest frame velocity of the star. We also observed a blue-shifted 
absorption component on 27 October superposed on the blue- 
shifted emission component, which is considered as a signature 
of the outflow. The full width at half maxima (FWHM) of the 
outflow increased with time, but its depth decreased (Table. [2}. 

3.6. Profile variations on 11, 28 and 29 December, 2008 

The Ha line profile on 11 December 2008 was of an in- 
verse P Cygni type with a narrow absorption centred at 
~30 km s^Fig.Wj). It was one of the few profiles where the red 
component almost matched the photospheric component. During 
28 December and 29 December 2008 the Ha line profiles were 
also of the inverse P Cygni type. We observed the emergence of 
an infall event on 28th December and followed its evolution for 
several hours and until the next night. The depth of the RAC 
was less than the depth of the central absorption component, 
which became red-shifted and wider compared to its counter 
part on 11 December. The nightly variations of the maximum 
velocity of the blue wings were between -250 to -300 km s _1 , 
and that of the red absorption wings was 300 km s~'and it did 
not vary much throughout the night. On 29 December 2008 
we observed both prominent infall and weak outflow events si- 
multaneously (Fig.gja) MHJD 830.100 to 830.280). The RAC 
that emerged on 28 December became deeper and its centre 
was separated from the central absorption component, as a re- 
sult of which the width of the central absorption component 
reduced (Fig. [9]). The maximum velocity at the blue wings did 
not change, but the average maximum velocity at the red wings 
shifted to -380 km s _1 on 29 December. The depth of the RAC 
varied throughout the observations on 28 and 29 December. We 
show the emergence and evolution of the RAC in Fig. [9] The ve- 
locity of the RAC was decelerated within a few hours from its 
starting point and then became constant at MHJD 830.21, where 
it remained until the end of the observation on 29 December 
(Table. [2]). All observed line profiles of 28 and 29 December 
showed a blue-shifted emission component whose strength in- 
creased after MHJD 829.1. We detected the outflow on top of 
the blue-shifted emission at MJHD 830.100 and could identify 
it clearly before MHJD 830.330. It gradually moved towards the 
central absorption features so the probable reason for its non- 
detection after MHJD 830.330 might be because the BAC came 
in the same line of sight of the central absorption. 
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Fig. 9. Selected Ha emission line profiles on 28 and 
29 December 2008 in the rest frame velocity of V35 1 Ori are 
overplotted. MHJD of the respective profiles are displayed using 
the same colour as of the profiles. 
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3. 7. Profile variations on 17 to 20 January, 2009 

The Ha emission line profiles observed during 17 to 
20 January 2009 are mostly dominated by both blue- and red- 
shifted emission components, i.e. mostly Type II and Type III 
profiles (Fig. [5}. The strength of both emission components var- 
ied during our observation, e.g. one of the component was al- 
ways weaker compared to the other component and thus the 
shapes of the profiles also changed from one type to another. The 

average maximum velocities at the wings were 250 km s _I 

and ~250 km s _1 in the blue and red sides respectively. We also 
detected a strong outflow, which was clearly detectable from 
17 to 19 January. The central narrow absorption component at 
~20 km s" 1 remained constant throughout the observations. We 
display the time evolution of the outflow in Fig. 10 The out- 



flow gradually became wider as it continuously moved towards 



Fig. 10. Selected Ha emission line profiles on 17 to 
20 January 2009 in the rest frame velocity of V351 Ori are 
overplotted. MHJD of the respective profiles are displayed us- 
ing the same colour as of the profiles. Vertical dashed lines of 
different colours represent the centres of the BACs estimated by 
eye on the respective MHJDs. Red-shifted emission increased on 
MHJD 852.070 after the disappearance of the outflow event. 



the central absorption. It is also worth to mention that the width 
of the central absorption was at its minimum among our ob- 
served spectra. The non-detection of the outflow on 20 January at 
MHJD 852.070 was most probably caused by the superposition 
of its line of sight with the central absorption. As a result, the 
central absorption became wider and asymmetric, e.g. it became 
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elongated toward blue side and the red-shifted emission reached 
its maximum strength compared to the other spectra observed 
during 17 to 19 January, 2009. 

3.8. Profile variations in March and April, 2009 

Most of the profiles we observed on several days in March 
and April were double-peaked profiles (Fig. |6|a)). On 2 and 
3 March 2009, there was a change in the slope of the blue-shifted 
emission. The red-shifted emission gradually gained strength 
and the profile evolved from type III-R to type II-R. The aver- 
age maximum velocity at the blue wing was 250 km s 

and -300 km s _1 at the red wing. We did not find any red- 
shifted emission on 29 March 2009 (MHJD 920.160). The pro- 
file showed a blue-shifted emission and the central absorption 
component. On the next day, we observed a strong red-shifted 
emission comparable to the blue one, and thus the profile was 
changed into a type II-R. The red-shifted emission was again 
considerably reduced on 3 1 March and we observed a type III- 
R profile. We detected a blue-shifted absorption component at 
~71 km s on 26 April 2009. The absorption component shifted 
towards the deep minima on 27 April. There was no signature of 
outflow in the spectra on 28 April. The red-shifted emission be- 
came gradually stronger from 26 to 28 April and became com- 
parable to the blue-shifted emission. These changes are similar 
to those observed on 28 and 29 December 2008 (Fig.Q, but the 
central absorption features were different, e.g. in April 2009 the 
central absorption component was board and asymmetric and ex- 
tended towards the red region compared to the narrow and nearly 
symmetric central absorption in December 2008. 

3.9. Hp, NaDl, and NaDl profiles 

We obtained simultaneous spectra containing the Hp, NaDl, and 
NaDl profiles along with Ha emission line profiles. Due to the 
poor signal-to-noise ratio, we could not analyse the individual 
Hp profiles. Instead, we plot the nightly average profiles of Hp 
in the rest frame velocity of V351 Ori and overplot the synthetic 
spectra in Fig. [TTJa). The shape of the Hp profiles also showed 
variation over the timescales of month. We did not find any emis- 
sion components in the Hp profiles. However, the prominent in- 
fall event occurred during 28 and 29 December, as can be easily 
seen also in the Hp line profiles. The NaDl and NaDl profiles 
displayed in Fig. [TT|b) consist of a broad and a narrow absorp- 
tion component. We did not find any elongated absorption in the 
broad line profiles. We measured the velocity of the narrow com- 
ponent and found that it was showing an almost constant velocity 
of ~ +22 km s _1 , i.e. red-shifted by ~1 1 km s _1 , compared to the 
star. 



4. Discussion 

HAeBe stars show photometric variabilities on timescales of 
days to years, [van den Ancker et aL] ( 1996) analysed the long- 
term optical photometric behaviour of V35 1 Ori and found that 
the star showed a large variation of >2 mag before 1985. |van| 
|den Ancker et al!] also reported that the visual magnitude of 
V351 Ori varied by 0.16 mag between 1985 to 1990, i.e. over 
timescale of years. During 1990 to 1995, the visual magnitude 
of the star became almost constant without any significant vari- 
ation. 

-0.46 mag) over 
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Fig. 11. (a) Nightly average Hp profiles of V351 Ori on several 
days, (b) Sample NaDl and NaDl profiles of V351 Ori on a 
few epochs. Velocities are in the rest frame of V351 Ori. 



variation within 0.05 mag. We obtained optical BVRI photom- 
etry of the star on 1 1 November, 2009 and obtained the mag- 
nitudes V-8.89, B-V -0.38, V - R -0.19, and V - I -0.41 
mag, which indicates that there are no significant changes in the 
long-term photometric behaviour of V351 Ori. However, it is 
also noteworthy that our photometric observation did not coin- 
cide with the spectroscopic observation. Unlike the variations 
in visual magnitudes, which reflect the variable stellar photo- 
sphere, near infrared (NIR) variability can account for changes 
in both the stellar photosphere and the circumstellar environ- 
Balona et al.| ( 200"2] > compared the NIR magnitudes of 



ments. 

V351 Ori and found thaTi/ - K excess of V351 Ori (which 
serves as a proxy for the inner disk materials) decreased con- 
tinuously over the period from 1982 to 2002. [yah den Ancker| 
|et al. ( |1996| l and Balona et al. (2002) suggested that the system- 
atic NIR variability of V35 1 Ori is caused by the clearing of dust 
around the star, which indicates a variable circumstellar environ- 
ment around V35 1 Ori. 

Our high-resolution spectroscopic observations revealed the 
variable circumstellar environment and presences of TACs in the 
weak Ha emission line of the Herbig Ae star V35 1 Ori. As men- 
tioned earlier, TACs in higher Balmer lines and in a few metallic 
lines have been detected in Herbig Ae stars such as BF Orionis, 
SV Cephei WW Vulpeculae, XY Persei, and UX Orionis ( [Mora] 
|et al-1 2002, 2004). All these Herbig Ae stars show prominent Ha 
emission profiles with equivalent widths r anging from -2.3 A 



(for UX Ori ) to -30.0 A (for WW Vul) ( |Manoj et al.l|2006} . 



Balona et al.| ( |2002| ) obtained the optical magnitudes (e.g. 
V-8.92, B-V -0.37, V - R -0.22 and V - 1 
timescales of minutes to weeks and did not find any significant 



These stars also show optical photometric and polarimetric vari- 
abilities over timescales of months (Oudmaijer et al. 200 1 ] > . 
XY Per, UX Ori, BF Ori, an d WW Vul also sho w NIR variability 
over timescales of months ( Eiro a~eT al. |200 1 } . The characteris- 
tics of V35 1 Ori do not match the trend of the Herbig Ae stars 
mentioned above, because V35 1 Ori shows nearly constant op- 



R. Choudhury et. al.: Variable circumstellar activity of V351 Ori 



9 



tical magnitudes over years and a low Ha emission equivalent 
width (< 1 A). Moreover, V35 1 Ori is also comparatively older 
(~6.5 Myr) than the mentioned group of Herbig Ae stars with 
ages ranging from 2.5 to 5.2 Myr (Manoj et al. 2006 Montesinos 
et al.|2009 1. These differences make V351 Ori an interesting tar- 



get to study the interaction of intermediate PMS stars (contract- 
ing towards the main sequence) with their circumstellar environ- 
ment. 

V351 Ori is known as an isolated HAeBe star without 
any associated nebulosity in its immediate vicinity. However, 
(Hernandez et al. 2005 1 suggested its association with the 
Ori OBlbc region located at a distance of 400 to 450 pc. The 
estimated ages of V35 1 Ori available in the literature vary from 
1 to 6.5 Myr ( |van den Ancker et al.fT996l |Ripepi et aLpOOl) . 
The main sources of uncertainty to estimate the age from the 
HR diagram is the distance to the star. The lower limit of the 
distance is ~210 pc (van Leeuwen 2001) and the upper limit 
is ~450 pc, which comes from the assumption of its associa- 
tion with the Ori OBlbc region (Hernandez et al. 2005). Ripepi 
et al. ( 2003) 1 analysed the pulsating behaviour of V35 1 Ori and 



suggested an intermediate value for the distance between the up- 
per and lower limit, which leads to an estimate of the age of 
~6.5 Myr. The Ha emission from HAeBes, which is a proxy 
for active accretion, are found to decrease significantly beyond 
3 Myr e.g. |Manoj et al.| ( |2006[ ). There are some HAeBe stars 
older than ~ 6 Myr with significant Ha emission (Manoj et al 
|2006 1 IMontesinos et al. ~~ 



2009). However, the proposed age of 



V351 Ori i.e. 6.5 Myr seems to be consistent with the common 
trend of weak Ha emission from stars with ages > 3 Myr. We 
also consider that the mass of V35 1 Ori is 2 M Q as estimated by 
Ripepi et al.|(|2003"i. 



4.1. Kinematics of TACs 

Blue-shifted absorption components (BACs) and red-shifted ab- 
sorption components (RACs) detected in Ha emission line pro- 
files of V351 Ori showed variation in the velocity centres, depths 
and widths of the profiles. We show the evolution of a RAC and 
BAC in Figs. [9] and [10] respectively. We tried to reproduce the 
observed Ha profiles by fitting a multiple number of Gaussian 
component but we also noticed that the fitted Gaussians were 
not unique. When we tried to fit some profiles, e.g. the ob- 
served profiles of 29 to 31 March, 2009, we found that we got 
a good fit when we fitted two separate Gaussians to the blue- 
shifted and red-shifted emission components rather than a sin- 
gle Gaussian for the whole emission component. The fitting of 
these complicated and variable profiles for all observed spec- 
tra in a consistent way demands a complete understanding of 
the geometry of line-forming regions and their distribution with 
respect to the central star, which is beyond the scope of this 
work. Accordingly we measured the velocity centres and depths 
of the TACs' profiles by fitting Gaussian profiles to the individ- 
ual TACs of the normalised spectra over the baseline of unity 
(obtained by the normalisation of the observed spectra by syn- 
thetic photospheric spectra). In this way, we estimated the veloc- 
ity centres and depths of the individual component confidently, 
however, we could not estimate the FWHM of the individual 
TACs, which requires the fitting of the entire profile. As a re- 
sult we discuss the kinematics of the observed TACs and not the 
physical conditions of the infalling and outgoing materials, be- 
cause that partly depends on the values of the FWHM as well. It 
is also worth to mention that the detection and analysis of TACs 
in other Balmer lines such as HB, Hy etc. provides additional 



information. Because of the poor signal-to-noise ratio at HB we 
could not detect the TACs as Ha profiles. However, good tem- 
poral coverage of the spectroscopic observations helped to make 
an unambiguous detection of the TACs in the Ha line profiles, 
and thus the kinematics of the TACs estimated from Ha pro- 
files only are also highly reliable. We list the velocity centres 



Table 2. Transient Absorption Components appeared in Ha line 
profiles of V35 1 Ori 
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and depths of the TACs from the normalised continuum level in 
Tab le [2] We plot the time evolution of the velocities of TACs in 
Fig 12 All detected TACs in the Ha line profiles showed decel- 
eration of the gas at various rates. The rate of deceleration of the 
outflow event on 27 October 2008 was ~2 m s~ 2 , for the infall 
event on 29 December 2008 it was ~1.7 m s~ 2 , for outflow event 



it was ~2.2 m s , and for the outflow event during the 17 to 
19 January 2009 it was fraction of m s~ 2 . Neither the magneto- 
spheric accretion model nor the wind model provide any detailed 
explanation about the deceleration of the TACs. A detailed cal- 
culation of a time-dependent model considering accretion and 
wind models and including the observed results is very impor- 
tant to understand the variable star-disk interactions. Broad red- 
shifted absorption wings at high velocities are supposed to orig- 
inate from the free fall of the surrounding material onto the star. 
The cool outflowing gas is supposed to originate farther away 
from the star than the accretion-dominated emission, but still 
within about -10 ^Edwards et al.||2003) . We found that the 
average maximum velocity of the red-shifted emission wings 
was ~250 to ~300 km s" 1 . We detected an enhanced red-shifted 
absorption wing extending up to ~380 km s~'on 29 December, 
2009. Assuming the free-fall velocity condition, we estimate the 
extent of the region involved in the accretion processes as ~5 R Q , 
which indicates that the changes are taking place within a radius 
of - 



-0.02 AU. |Monnier & Millan-Gabet| p002 ) studied the dust 
sublimation radius of the PMS stars with a stellar radiation field. 
They discussed the relation between the size of the inner dusty 
disks and the luminosity of the central stars in the context of 
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computed using the relation 
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Fig. 12. Time evolution of the TACs as they appeared in the 
Ha line profiles. Triangles represent the outflow event on 
27 October 2008, squares and circles represent the outflow and 
infall events respectively on 29 December 2008, and diamonds 
represent the outflow event during 17 to 19 January 2009. 



different disk models including the paradigm of both optically 
thick and thin inner disks. Considering the uncertainties such as 
the distance and thus the luminosity, and the model for optically 
thin/thick disks, the typical dust sublimation radius of V35 1 Ori 
would be in the range of ~ 0.07 to .2 AU. A comparison be- 
tween the region of the gas accretion and the dust sublimation 
radius hint towards the presence of gaseous material in the in- 
ner disk cavity of V351 Ori. Further investigation of the inner 
disk properties of PMS stars such as V35 1 Ori will be helpful 
to understand the disk-dissipation processes and the associated 
timescales. 



4.2. Dust mass estimates using IRAS fluxes 

The detection of red-shifted absorption wings that extend be- 
yond 300 km s _1 supports a scenario of the accretion of sur- 
rounding materials onto the star. However, as V351 Ori also 
shows weak Ha emission, so it will be interesting to investigate 
whether the left-over gas in the circumstellar environment of the 
star is sufficient to continue the accretion. We estimate the dust 
mass and thus the gas mass in the circumstellar disk of V351 Ori 
using the Infrared Astronomical Satellite (IRAS) fluxes. 

Dust grains present in the circumstellar environment absorb 
and scatter varying fractions of the incident photons (A < 1/j.m) 
and re-emit part of the absorbed energy in the infrared band. Far 
infrared emission properties of the dust grains can be used to 
derive estimates on the amount of dust of the disks of PMS stars, 
assuming the emission at these wavelengths is optically thin. As 
suggested by Hildebrand] ( |1983| l, the disk dust mass M d can be 



M d 



4tt ap d d 2 F v (T d ) 
3 Q v (T d ,a) nB v (T d ) ' 



where a, p d and d are the grain radius, specific grain mass 
density, and distance to the source, respectively. F v , Q v (T d ,a), 
and B v (T d ) are the observed flux density, grain emissivity and 
the Planck function of temperature Td at frequency v. |van den| 
|Ancker et al. ( |1996[ ) analysed the spectral energy distribution 
(SED) of V35 1 Ori and estimated that the dust temperature cor- 
responds to the far infrared (FIR) IRAS fluxes as 80 K. The 
FIR SED of V35 1 Ori peaks between 25 and 60 fim and drops 
gradually between 60 and 100 fan, which suggests that most 
of the emission arises beyond 25 yum. In that case the calcula- 
tion of dust grain temperature using 60 and 100 pm flux den- 
sities ( |G. Helou & D. W. Wa lker 1988) would be a bet ter es- 

i s \"- 4 J 

timate. Using the following relation, Tj = 49 [-^J (Young 
et al.||1989] l, we estimated the dust grain temperature as 53 K. 



and 



= 0.024 and 



We adopted p d — 3 gem 3 , a = 0.1 fim 

0.04 gem' 2 at 60 and 100 jim respectively (Hildebrand 1983). 



2&_ 



The estimated dust masses at 60 and 100 /urn are 9.5x10^ and 
9.2xl0~ 5 M Q respectively at an adopted distance of 450 pc and 
2.07xl0~ 5 , and 2.01xl0~ 5 M at 210 pc. Assuming a gas-to- 
dust ratio of 100, the estimated disk mass (dust+gas) would be 
in the range of 0.002 to 0.009 M Q . The estimated values of the 
disk mass of V35 1 Ori are lower compared to the disk mass of 
actively accreting Herbig Ae type stars such as AB Aur (0.01 
M Q ), H P 163296 (0.028 Af Q ) deri ved from millimetre obser- 
vations (Mannings & Sargen t]! 997} . This is consistent with the 
relatively evolved circumstellar environment of V351 Ori. 



4.3. Accretion and wind models 

The Ha emission line serves as one of the most empirical sig- 
natures of accretion of the circumstellar material onto the cen- 
tral PMS stars. Magnetospheric accretion models are quite suc- 
cessful in explaining the observed spectral characteristics such 
as broad asymmetric line profiles etc. of accreting classical T 
Tauri stars (CTTS) i.e. the lower mass counterparts of HAeBe 
Partmann et al.|1994||Muzerolle et al.|1998||2001| >. According 
to magnetospheric accretion models, the stellar magnetic field 
disrupts the circumstellar disk at several stellar radii and the ac- 
creting material falls onto the stellar surface along the magnetic 
field lines Koenigl ( 1991 ). The results of magnetospheric accre- 
tion models suggest that the hydrogen emission lines form in 
the infall zone, so that the blue-shifted asymmetric emission-line 
profiles arise due to partial obscuration of the flow by the inner 
part of the accretion disk and the red-shifted absorption profiles 
result from infalling material at near free-fall velocities onto the 
stellar surface ( |Muzerolle et al.||1998"] l. jMuzerolle et al.| ( |2004) > 
also demonstrated the magnetically channelled disk accretion in 
the Herbig Ae star UX Ori and partly reproduced the double- 



peak Ha profile of UX Ori. Several studies e.g. Alencar et al. 
( |2005| l; [Kurosawa et al.| ( [2006) l found that disk wind also play an 
important role in the formation of double-peak profiles in CTTS. 

V351 Ori showed varieties of Ha line profiles throughout 
our observation. The shape of the profiles contains informa- 
tion about the geometrical orientation and physical processes 
involved in the interaction of the star with its circumstellar envi- 
ronment. The common trends in all spectra are the slightly red- 
shifted central absorption and the elongated blue-shifted emis- 
sion components. The width of the central absorption changed 
within a day. Though the strength of the blue-shifted emission 
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component varied with time, but there was no significant vari- 
ation in the maximum velocity at the wings. On the contrary, 
the red-shifted component showed drastic variations in shape, 
strength, and velocity. The strength of the red-shifted emission, 
if present, was less than the blue-shifted emission, which can 
be attributed to the selective absorption of the red-shifted flux 
due to the asymmetric distribution of the star-disk system with 
respect to our line of sight. 



Kurosawa et al.| (|2006)> investigated the formation of Ha 



emission line profiles around CTTS combining the magneto- 
spheric accretion and disk-wind models (hybrid model). Though 
the Ha emission line profiles are calculated with typical param- 
eters of CTTS, but it has also been suggested that the model 
can be used for other PMS stars such as HAeBe stars and for 
Brown dwarfs (BDs). In order to understand the relative impor- 
tance of the accretion and wind on the formation of the line pro- 



file Kurosawa et al. [ calculated the Ha line profiles for three cases 
e.g. the magnetospheric accretion, the disk wind, and a combi- 
nation of magnetospheric accretion and disk wind models. With 
the use of this hybrid model they were able to reproduce the 
variety of the Ha emission line profiles reported by Reipurth 
|et al.|([T996]l.|Kurosawa et aTT also noticed the degeneracy in the 



line profiles that emerged from different models mostly from the 
disk wind model and the hybrid model. For example, the dou- 
ble peak Ha emission line profiles can be reproduced using the 
disk-wind model only. It is somewhat difficult to distinguish the 
contribution of magnetospheric accretion and disk winds in the 



emission line profiles . Alencar et al. ( 2005 1 reported the line 



profile variability of CTTS RW Aur and found significant varia- 
tions in the intensity of the blue- and red-shifted emission com- 
ponents similar to V351 Ori. Alencar et al. tried to reproduce 
the double-peaked Ha emission line profiles with the magneto- 
spheric accretion model and the disk wind model. The Ha, Ha, 
and NaDl profiles calculated using the disk wind models pro- 
vided a better fit to the observed spectra. These studies reveal 
the importance of disk wind in the emission line profiles of the 
PMS stars. While the magnetospheric accretion model produces 
double-peak profile under certain conditions such as accretion 
rates etc., the wind model produces double-peak profiles in usual 
conditions (e.g. |Alencar et al.| ( |2005[ ); |Kurosawa et"aL] ([2006)). 
So it is possible that the disk wind contributes significantly in 
the Ha emission of V35 1 Ori, which showed a double-peak pro- 
file most of the time during our observation. However, we did not 
detect [O I] /16300 A emission in our observed spectra, which ac- 
cording to sever al authors forms in the disk wind (Hartigan et al. 
Il995l. However, 



Acke et al. 



(2005 ) proposed that [O I] /16300 A 



emission arises in the surface layers of the protoplanetary disks 
surrounding HAeBe stars due to the the photodissociation of OH 
molecules. The maximum equivalent width of the Ha emission 
in our observed spectra is < 1 A, which is less than the typical 
minimum equivalent width of the Ha emission, of ~4 A of the 



observed samples of HAeBe stars studied by |Acke et a l. (2005). 
Considering the discrepancy regrading a thermal or non-thermal 
origin of the [O I] /16300 A emission and the low Ha emission, 
the absence of the line may not contradict with the presence of 
a disk wind. However, considering the rapid variation in the Ha 
emission on timescale of a day, which can be attribute to the 
variation of wind acceleration rate (Kurosawa et al. 2006}, we 
suggest that the disk wind may be responsible for most of the 
Ha emission of V351 Ori. However, the degeneracy of accre- 
tion and disk-wind models cannot be explored completely with 
our observed spectra. 



We try to interpret the behaviour of V35 1 Ori in the context 
of hybrid model, i.e. by considering the effect of magnetospheric 
accretion as well as the disk-wind. The central absorption that 
was present in most of the spectra seemed to arise because of 
the geometrical configuration of the star-disk system which is 
viewed viewed nearly edge-on. The changes in the shapes of the 
central absorption component, for most of the epochs, seemed 
to arise from the superposition of infall or outflow events in 
the same line of sight. On 27 October 2008, the Ha emission 
line profile may have been powered by a disk wind with the si- 
multaneous presence of infall and outflow. The infall was quite 
prominent and it obscured the red-shifted emission component. 
Before 30 October 2008, the infall event ended and the central 
absorption returned to the normal shape with a nearly symmet- 
ric double-peak profile. On 28 and 29 December 2008, the Ha 
emission line profiles were dominated by the infall event. During 
17 to 20 January 2009, the Ha emission line profiles were dom- 
inated by the disk wind but the blue-shifted emission was par- 
tially obscured by the prominent outflow event. The Ha emis- 
sion line profiles observed during March and April 2009 also 
dominated by the disk wind. The Ha emission from V351 Ori 
seems to originate in disk wind and at the same time in sporadic 
infall of matter moderated by a magnetic field, and the outflow 
of cool gas occurs quite frequently. 



5. Conclusion 

Spectral analysis of 45 high-resolution spectra reveals that the 
line profile variability of V35 1 Ori is ubiquitous. The line profile 
variation originate due to the disk wind, rate of wind accelera- 
tion, episodic accretion, and cool outflowing gaseous material. 
These signatures support the complex and dynamical interaction 
of the circumstellar material with V351 Ori. The circumstellar 
environment of the Herbig Ae star V35 1 Ori consists of an inner 
cavity with gaseous material, a dusty disk and disk wind. The 
interaction processes are highly time-dependent from one hour 
for the changes in TACs to one day and longer timescales of 
the order of months for the changes in the overall characteristics 
of the broad Ha profiles. Owing to the extensive time coverage 
of the observation down to timescales of hours we could fol- 
low the time evolution of the individual TACs unambiguously. 
The long term monitoring of the star at timescales of months 
confirmed the deviation of the Ha emission line profiles from 
inverse P Cygni types. Accretion and outflow events are not con- 
tinuous in nature, and the profiles of the TACs seem to be cre- 
ated by obscuration due to gaseous blobs, which indicate that 
the circumstellar environment of V351 Ori is inhomogeneous 
and clumpy in nature. Detectable changes in the TACs showed a 
deceleration of the order of a fraction to a few m s -2 . Variation 
in the blue- and red-shifted emission components also occurred 
frequently on the shortest timescale of hours, which suggests the 
rate of the wind acceleration is also time-dependent. Finally, the 
presence of emission in the absence of veiling emphasises the 
importance of the disk wind in V35 1 Ori. Models containing dy- 
namic and non-axysymmetric magnetospheric accretion with the 
disk wind may be able to provide a satisfactory explanation of 
the complex circumstellar activities of PMS stars. 
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